A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

06(08.03.4;08.05.2;08.09.2 4U 0115+63;08.02.1;08.14.1;13.25.5) 



ASTRONOMY 

AND 
ASTmPHYSICS 



The Be/X-ray transient 4U 0115+63/V635 Cas 



I. A consistent model 



O 
O 

o 

(N 
> 

o 

(N 



> 
O 



o , 
o ■ 

0,; 

6 ' 
a: 



Ignacio Negueruela 1,2 * and Atsuo T. Okazaki 3 

1 SAX SDC, ASI, c/o Nuova Telespazio, via Corcolle 19, 100131 Rome, Italy 

2 Astrophysics Research Institute, Liverpool John Moores University, Byrom St.. 

3 Faculty of Engineering, Hokkai-Gakuen University, Toyohira-ku, Sapporo 062-8605, Japan 

Received / Accepted 



Liverpool, L3 3AF, UK 



Abstract. We present photometry and high SNR spec- 
troscopy in the classification region of V635 Cas, the 
optical counterpart to the transient X-ray pulsator 
4U 0115+63, taken at a time when the circumstellar enve- 
lope had disappeared. V635 Cas is classified as a B0.2Ve 
star at a distance of 7 — 8 kpc. We use the physical pa- 
rameters derived from these observations and the orbit 
derived from X-ray observations to elaborate a model of 
the system based on the theory of decretion discs around 
Be stars. We show that the disc surrounding the Be star 
must be truncated by the tidal/resonant interaction with 
the neutron star and cannot be in a steady state. This 
explains many of the observed properties of 4U 0115+63. 
In particular, because of this effect, under normal circum- 
stances, the neutron star cannot accrete from the disc, 
which explains the lack of regular Type I outbursts from 
the source. 
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1. Introduction 

The hard X-ray transient 4U 0115+63 (X 0115+634) is one 
of the best studied Be/X-ray binary systems (see Cam- 
pana 1996; Negueruela et al. 1997, henceforth N97). More 
than 50 of these systems, in which a neutron star orbits 
a Be star in a moderately eccentric orbit, are known (see 
Negueruela 1998; Bildsten et al. 1997). The Be star is sur- 
rounded by a disc of relatively cool material, presumably 
ejected from the star due to causes unknown, but gener- 
ally believed to be associated with fast rotation, magnetic 
fields and/or non-radial pulsations (Slettebak 1988). The 
presence of the disc gives rise to emission lines in the op- 
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tical and infrared spectral regions and an excess in the 
infrared continuum radiation. 

The hard X-ray emission is due to the accretion of 
circumstellar material on to the neutron star companion. 
Due to their different geometries and the varying physi- 
cal conditions in the circumstellar disc, Be/X-ray binaries 
can present very different states of X-ray activity (Stella 
et al. 1986). In quiescence, they display persistent low- 
luminosity (L x <S 10 36 erg s _1 ) X-ray emission or no de- 
tectable emission at all. Occasionally, they show series of 
periodical (Type I) X-ray outbursts (L x w 10 36 — 10 37 erg 
s^ 1 ), separated by the orbital period of the neutron star. 
More rarely, they undergo giant (Type II) X-ray outbursts 
(L x <; 10 37 erg s _1 ), which do not clearly correlate with 
the orbital motion. Some systems only display persistent 
emission, but most of them show outbursts and are termed 
Be/X-ray transients. 

The transient 4U 0115+63 was first reported in the 
Uhuru satellite survey (Giacconi et al. 1972; Forman et al. 
1978), though a search of the Vela 5B data base revealed 
that the source had already been observed by this satellite 
since 1969 (Whitlock et al. 1989). Precise positional deter- 
minations by the SAS 3, Ariel V and HEAO-1 satellites 
(Cominsky et al. 1978; Johnston et al. 1978) were used to 
identify the system with a heavily reddened Be star with 
a visual magnitude V » 15.5 (Johns et al. 1978; Hutch- 
ings & Crampton 1981), which was subsequently named 
V635 Cas (Khopolov et al. 1981). Rappaport et al. (1978) 
used SAS 3 timing observations to derive the orbital pa- 
rameters of the binary system, which consists of a fast- 
rotating (P s = 3.6 s) neutron star in a relatively close 
(-Porb = 24.3 d) and eccentric (e = 0.34) orbit around the 
Be star (see also Tamura et al. 1992). 

Due to the fast rotation of the neutron star, centrifugal 
inhibition of accretion prevents the onset of X-ray emis- 
sion unless the ram pressure of accreted material reaches 
a relatively high value (Stella et al. 1986; N97). The sys- 
tem had only be known to display Type II activity until a 
short series of Type I outbursts was detected by BATSE 
and RXTE in 1996 (Bildsten et al. 1997; Negueruela et 
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Fig. 1. The spectrum of V635 Cas in the classification region. Two exposures taken on November 14, 1997, with ISIS 
on the WHT equipped with the R1200B grating and the EEV10 camera have been combined for this figure. The 
comparison spectrum is that of the B0.2V standard r Sco. Both spectra have been divided by a spline fit to the 
continuum for normalisation and smoothed with a a = 0.8A Gaussian function for display. 



al. 1998). The giant outbursts are associated with large 
amplitude brightenings of the optical and infrared magni- 
tudes of the counterpart (N97). 

This is the first of two papers dedicated to providing a 
coherent picture of this system and understanding the im- 
plications of its unusual behaviour for the general class of 
Be/X-ray transients. Here we derive the astrophysical pa- 
rameters of 4U 0115+63 and build a model for circumstel- 
lar disc around the Be star that allows us to understand 
the usual quiescent state of the system. In the second pa- 
per (Negueruela et al. 2000; henceforth Paper II), we will 
analyse the temporal evolution of this disc and investigate 
how its behaviour is connected with the X-ray activity of 
the source. 

2. Observations 

We present optical photometry and optical spectroscopy of 
V635 Cas, obtained, for the first time, when the emission 
component was almost completely absent. 

2.1. Optical spectroscopy 

Observations of the source were taken on November 14, 
1997, using the Intermediate Dispersion Spectroscopic and 
Imaging System (ISIS) on the 4.2-m William Herschel 
Telescope (WHT), located at the Observatorio del Roque 
de los Muchachos, La Palma, Spain. The blue arm was 
equipped with the R1200B grating and the EEV10 CCD, 



which gives a nominal dispersion of ~ 0.22 A/pixel over 
~ 900 A. The resolution at ~ A4600A, estimated from the 
FWHM of arc lines, is ~ 0.7 A. Two exposures were taken, 
centred on A4250A and A4650A. A combined spectrum is 
shown in Fig. |]. 

The red arm was equipped with the R1200R grating 
and the Tek5 CCD, which gives a nominal dispersion of 
~ 0.4 A/pixel at Ha (the resolution is ~ 0.8A at Ha). A 
similar observation has been attempted on July 19, 1997, 
but due to bad weather conditions, the blue spectrum is 
too noisy to be of any use. An Ha spectrum, though, was 
obtained using the same setting as in the November obser- 
vations. The observations are displayed in Fig. 0. All the 
data have been reduced using the Starlink software pack- 
ages CCDPACK (Draper 1998) and FIGARO (Shortridge et 
al. 1997) and analysed using FIGARO and DIPSO (Howarth 
et al. 1997). 

2.2. Optical photometry 

From the observations listed in Paper II, we have taken the 
dataset showing the faintest and bluest magnitudes, which 
corresponds to 1998 January 7. The values measured are 
U = 17.05 ± 0.10, B = 16.92 ± 0.07, V = 15.53 ± 0.04, 
R = 14.55 ± 0.03 and I = 13.49 ± 0.03. As discussed 
in Paper II, these observations coincided with a disc-less 
phase of V635 Cas and should represent values close to the 
intrinsic magnitudes of the star. We note that the value 
for / is compatible with the faintest point measured by 
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Fig. 2. High resolution spectra of V635 Cas showing the 
almost completely disc-less state observed in 1997. In the 
spectrum from July 19 (dashed line) only some residual 
emission is visible in the He I A6678 A line. In the Novem- 
ber 14 spectrum, taken simultaneously with the blue spec- 
troscopy, the disc is beginning to reform and some emis- 
sion is present now in Ha. The peak separation corre- 
sponds to Aup ea k ~ 600 km s -1 , well above what would be 
seen in the case of angular momentum conservation. The 
strong narrow features are interstellar diffuse bands. The 
spectra have been divided by a spline fit to the continuum 
for normalisation. 



Mendelson & Mazeh (1991) in the five years that their 
photometry covers. 

3. Results 

3.1. Spectral classification 

Unger et al. (1998) have presented intermediate resolu- 
tion spectroscopy of V635 Cas and classified the source 
as 09.5V. However, their spectrum had a very low signal 
to noise ratio and showed several spurious features. The 
spectrum shown in Fig. [j] has a much higher SNR and 
better resolution. At the time it was taken, Ha was start- 
ing to develop emission peaks after several months during 
which it had been observed in absorption (see Paper II). 
All the other Balmer lines are still in absorption, but they 
all present a similar asymmetric shape, which must be 
due to new emission components starting to develop. Sur- 
prisingly, many He I lines show emission components (AA 
4009, 4026, 4144 and 4713 A), as can be easily seen by 
comparing the spectrum with that of the B0.2V standard 
r Sco (Fig.g). 

The broad shallow lines are typical of a main-sequence 
early star. Hen A4686 A can be clearly seen, though it is 
weak, indicating a spectrum earlier than BO. 7. The weak- 
ness of the Sim AA 4552, 4568, 4575 A triplet (if it is at 
all present) supports the early spectral type and indicates 
a main sequence star. No obvious On lines are present, 
though the Cm AA 4072, 4650 A lines could be blends. If 
the weak feature at A454lA is real, then V635 Cas could 



be as early as B0V, but not earlier - the Si IV lines, which 
must be hidden on the very broad wings of HJ, should be 
clearly visible for an earlier (or higher luminosity) star. 
Therefore the star is constrained to belong to the BO- 
BO. 5 range. For simplicity, we will adopt the intermediate 
B0.2Ve spectral type. 

3.2. Distance 

As discussed in Paper II, during the second half of 1997 the 
source showed Ha and all other lines in the red in absorp- 
tion. On January 7th, 1998, we measured the faintest and 
bluest magnitudes in our dataset. Here we investigate the 
possibility that these photometric magnitudes are close to 
the actual apparent magnitudes of the star without any 
contribution from a disc, implying that the measured red- 
dening is (almost) purely interstellar. 

The intrinsic colour of a B0V star is (B — V)o = —0.26 
(Wegner 1994) and therefore the measured (B — V) = 1.39 
implies a reddening E(B — V) — 1.65 ± 0.08. Using the 
relationship E{U -B) = E(B-V) [0.69 + 0.04 x E(B-V)] 
from Fitzpatrick (1999), we find E(U - B) = 1.25. The 
measured (U - B) = 0.13 implies (U - B) = -1.12 ± 
0.12, perfectly compatible with the expected (U — B) — 
— 1.08. Further assuming the standard extinction curve 
for R — 3.1 after Fitzpatrick (1999), we find values for 
(V - R) = -0.31 and (V - 7) = -0.60. These values 
are, within the errors, close to the intrinsic colours of a 
B0V star, though slightly too blue. Wegner (1994) gives 
(V - R) = -0.14 and (V - I) = -0.37. Such an effect 
is surprising, since any residual emission is expected to 
contribute more strongly at longer wavelengths. 

We have considered the possibility that the extinction 
law to V635 Cas is different from the standard, since the 
line of sight to V635 Cas could intersect the outer regions 
of the OB association Cas OB7. This association is about 
3° in extension, centred at ~ (/ = 123, b = +1) and the 
de-reddened DAI = 11.3 implies a distance of 1.8 kpc 
(Garmany & Stencel 1992). We have dereddened the pho- 
tometry using the extinction curves of Fitzpatrick (1999) 
for different values of R. We find that values of R different 
from the standard value by more than ~ 0.2 provide dered- 
dened energy distributions that are farther away from the 
theoretical values than that dereddened with R = 3.1. On 
the other hand, the derived colours become very close to 
those expected when a value of EiB — V) close to the 
lower limit allowed by the errors is considered. 

We have also measured the EW of the interstellar dif- 
fuse band at A6613A in several Ha spectra of V635 Cas 
(N97; Paper II). We have obtained an average value EW — 
480+100 mA. According to the relation by Herbig (1975), 
the implied colour excess is E(B — V) — 1.9 ± 0.4, con- 
sistent with the photometric value. Likewise, the EW of 
the interstellar diffuse band at A4430A is 4.4 ± 0.4, which 
results (Herbig 1975) in E(B - V) = 2.0 ± 0.2. We note 
that the Nai H and K lines cannot be used for distance 
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Fig. 3. A Kurucz model for T e g = 26000K and logg — 4.0 (the adopted parameters for V635 Cas) has been normalised 
to the .B-band magnitude and is compared to the measured energy distribution (points have been set at the effective 
wavelengths of the different bands) dereddened with R — 3.1 and E(B — V) = 1.60 (a value close to the lower limit 
allowed by the error bars in our photometry). 



calibration, since the values measured are well above the 
saturation value for the relationship with reddening (Mu- 
nari & Zwitter 1997). 

Even though the interstellar diffuse bands seem to sup- 
port a value for E(B — V) close to the upper limit allowed 
by the error bars, the photometric evidence supports a 
value close to the lower limit, as can be seen in Fig. [|. 
The choice of a lower E(B — V) results in a higher dis- 
tance for V635 Cas. 

Assuming an absolute magnitude My = —4.1 after 
Vaccaet al. (1996), and taking E(B-V) = 1.65±0.08, we 
derive a distance to V635 Cas of 8±lkpc. This could be an 
overestimate if some residual emission was still contribut- 
ing to the intrinsic luminosity of V635 Cas. Moreover, the 
calibration of Vacca et al. (1996) indicates that early-B 
stars are hotter than previously thought and therefore the 
intrinsic (B — V)o used is probably too low. We note that 
distances ~ 7 kpc are obtained for the Be/X-ray transient 
V0332+53 (Negueruela et al. 1999) and for the standard 
High Mass X-ray Binary 2S 0114+65 (Reig et al. 1996), 
which is very close in the sky to 4U 0115+63. We consider 
very likely that the distance to 4U 0115+63 is also rj 7 
kpc and that all these objects belong to an outer arm of 
the Galaxy, situated beyond the Perseus arm at a distance 
<; 4 kpc (Kimeswenger & Weinberger 1989). At this dis- 
tance, the maximum luminosity of the source, measured 



by Vela 5B during a Type II outburst in 1974 (Whitlock 
et al. 1989) is close to 10 38 erg s _1 , i.e., close to the 
Eddington luminosity for a neutron star. 

4. A model for 4U 0115+63 

4-.1. Stellar and orbital parameters 

Given the B0.2Ve spectral type of the optical component, 
we expect it to have a mass Af* ~ 19 Mq (Vacca et al. 
1996). In this case, the mass function f(M) — 5.0 Mq im- 
plies an inclination angle i = 42° for a standard neutron 
star mass M x = 1.4 Mq. It is noteworthy that making the 
companion star moderately undermassive will not change 
strongly the inclination angle. For example, if M» = 
12 Mq, the inclination angle is i = 54°. For any reasonable 
value of M* , 40° < i < 60° or 0.7 < sini < 0.8. The depen- 
dence on M x is still smaller. Therefore, in what follows we 
will adopt as a model M* = 18 Mq and M x = 1.4 Mq, im- 
plying i = 43°. In this case, the measured a x sini implies 
a binary separation a = 6.6 x 10 10 m = 95 i?Q. Assuming 
i?* = 8 Rq (Vacca et al. 1996), this results in periastron 
and apastron distances a por = 8 i?» and a ap = 16 i?*. 

We have tried to estimate the inclination of the equa- 
torial plane of V635 Cas by using Buscombe's (1969) ap- 
proximation and the fits to the FWHM of He I AA 4388, 
4471 A of Steele et al. (1999). This is only a coarse approx- 
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imation, since weak emission components are present on 
the wings of the lines and, given the relatively low SNR of 
the blue spectrum and the broadness of the Balmer lines, 
several weak lines could be contaminating the line shapes. 
In spite of this, all the measurements obtained fall in the 
range vsini 240 — 340 km s _1 with most of them con- 
centrated around vsini = 290 km s _1 and the Balmer 
lines consistently giving higher values than the Hei lines 
(except Ha). 

Therefore, allowing for the effect of emission contami- 
nation, we will adopt as an estimate vsini sa 300 km s _1 , 
taking into account that the value is not likely to be much 
smaller, but could be larger. We note that, if we assume 
coplanarity between the equatorial and orbital plane, this 
value implies a rotational velocity v — 440 km s _1 , which 
is compatible with the assumption that Be stars rotate at 
& 80% of their break-up velocity. Moreover, the double- 
peaked emission lines which are seen soon after the cir- 
cumstellar disc starts reforming are typical of a moderate 
inclination (i w 40° — 60°), confirming that the value of 
sini cannot be very distant. 

Therefore, we find no compelling reason to argue for 
an inclination between the orbital and equatorial plane. 
As discussed in Paper II, this does not imply that the 
circumstellar disc must always remain coplanar. 

4-2. Viscous disc around V635 Cas 

To date, discs around seven Be stars have been spatially 
resolved with optical interferometers (Quirrenbach et al. 
1997). The size of the Ha-emitting region Ru a of these 
stars ranges from 3 to 12 stellar radii. The separation of 
the two peaks of Ha, Aw pca k, is significantly wider than 
2vsini (i?Ha/-R*) -1 expected in an angular-momentum 
conserving disc, but is consistent with a disc structure 
in which the rotation velocity is near Keplerian. 

Although there is no widely-accepted model to form 
near-Keplerian discs around Be stars, a viscous decretion 
disc model proposed by Lee et al. (1991) seems promis- 
ing (Porter 1999; see also Okazaki 2000a). In this model, 
the matter supplied from the equatorial surface of the star 
drifts outward because of the viscous effect and forms the 
disc. Basic equations for viscous decretion discs are the 
same as those for viscous accretion discs, except that the 
sign of M (mass decretion/accretion rate) is opposite. The 
viscous decretion disc model thus predicts a geometrically 
thin and near-Keplerian disc around a Be star. The out- 
flow in the viscous decretion disc is highly subsonic near 
the star (Okazaki 1997, 2000b), which is consistent with 
the observed upper limit on the radial velocity (3kms _1 ) 
in discs around Be stars (Hanuschik 2000). 

Since the mass of the neutron star is much smaller 
than that of the Be star, its presence will hardly affect 
the mass loss process from the Be star and will have a 
minor effect on the formation and structure of the Be disc 
except in the region near the critical Roche potential. We 



Table 1. Model parameters for V635 Cas/4U 0115+63. 



Parameter 


Notation 


Value 


Binary orbit 






Period 


Porb 


24.3 d 


Inclination 


i 


43° 


Semi-major axis 


a 


95 R®(= 12 R,) 


Eccentricity 


e 


0.34 


Neutron star 






Mass 


M x 


1.4 Mq 


Be star 




Mass 


M t 


18.0 Mq 


Radius 


R, 


8.0 R® 


Effective temperature 


Peff 


26,000 K 


Be disc 






Temperature 


T d 


20,800 K 


Scale-height at r = R* 


H(R,)/R, 


2.6 ■ 10~ 2 



note that the overall characteristics of the discs in Be/X- 
ray binaries, as derived from the observation of emission 
lines and their evolution are not fundamentally different 
from those of isolated Be stars (see Negueruela et al. 1998, 
where a discussion of the similarities and main differences 
is presented; and see also Section ||). 

Therefore, as a first approximation, we will assume in 
this paper that the viscous decretion disc model, which 
has been developed as a disc model for isolated Be stars, 
is applicable to Be discs in Be/X-ray binaries. That is, 
we will assume that the disc surrounding V635 Cas is 
near-Keplerian when unperturbed. For simplicity, we also 
assume that the disc is isothermal with a temperature 
of 0.8T e ff, where the effective temperature of the star is 
T cff = 26000 K. Table @ summarizes the model parame- 
ters adopted for 4U 0115+63. 

4-3. Tidal truncation of the viscous disc 

Artymowicz & Lubow (1994) investigated the 
tidal/resonant truncation of the accretion discs in 
eccentric binary systems and found that the disc size 
becomes smaller in a system with larger eccentricity. 
Following their formulation of the companion-disc inter- 
action, we study below the tidal truncation of the Be-star 
disc in 4U 0115+63. 

In order to evaluate the tidal/resonant torque exerted 
by the neutron star on the Be-star disc, we consider the 
binary potential in a coordinate system (r, 6, z) in which 
the origin is attached to the Be star primary: 

s>( o ) GM * ^ 

l ' ,Z> ~ r [r 2 + 7-2 - 2rr 2 cos(0 - f)} 1 / 2 

+^cos(0-/), (1) 

where M* and M x are masses of the Be and neutron stars, 
respectively, r 2 is the distance of the neutron star from 
the primary, and / is the true anomaly of the neutron 
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star. The third term in the right hand side of Eq.([j]) is the 
indirect potential arising from the fact that the coordinate 
origin is at the primary. 

We expand the potential by a double series as 



$(r, ) z)=y~] (f) m i exp[j(m6» - Kl B t)], 



(2) 



m,l 



where m and I are the azimuthal and time-harmonic num- 
bers, respectively, and il B = [G(M* + M x )/a 3 ] 1 / 2 is the 
mean motion of the binary with semimajor axis a. The 
pattern speed of each potential component is given by 
O p = (l/m)£l B . 

Inverting Eq.(^) and denoting the angle 6 — f as ip, we 
have 



GM f 2 a 

= - i -o / d(tt B t)—cos(mf -lfl B t) 

la r 2 

cos mtp 



a 



X / d ^(l + /32-2/3cos^)V2 

-^1 f d(n B t)—0cos(mf - m B t)\ , (3) 
I" Jo r 2 J 

where (3 — r jr-i and S m \ is the Kronecker delta function. 

For each potential component <f> m i , there can be outer 
and inner Lindblad resonances at radii where £l p = 
f2 ± k/to and a corotation resonance at the radius where 
flp = O. Here, k is the epicyclic frequency, and here and 
hereafter the upper and lower signs correspond to the 
outer Lindblad resonance (OLR) and inner Lindblad res- 
onance (ILR), respectively. The radii of these resonances 
are given by 



and 



to ± 1 
I 



2/3 



(i + s) 



-1/3, 



/TO\ 2 / 3 . 

rcH=( T ) (1 



-1/3 



a, 



(4) 



(5) 



where g = M x /M*. 

For near-Keplerian discs, where f2 ~ k ~ 
(GM,/r 3 ) 1/2 , Goldreich & Tremaine's (1979; 1980) stan- 
dard formula for torques T m i at the outer Lindblad res- 
onance (OLR) and inner Lindblad resonance (ILR) is re- 
duced to 



± 



to(to ± l)n 2 a(X + 2m 
3/ 2 fil| 



(6) 



where a is the surface density of the disc at the resonance 
radius and A = (dlr±<^ m //<ilnr)LR- Similarly, the torque 
at the corotation resonance (CR) is written as 

_ 2m 3 ir 2 <jtf nl 

Note that angular momentum is removed from the disc at 
the ILRs, whereas it is added to the disc at the OLRs and 
CRs. 



In near Keplerian discs, the viscous torque formula de- 
rived by Lin & Papaloizou (1986) is written as 



I ' H 

T v i, = 3TraGM*crr ( — 
r 



(8) 



[see also Artymowicz & Lubow (1994)], where a is the 
Shakura-Sunyaev viscosity parameter and H is the verti- 
cal scale-height of the disc. In our isothermal disc model 
(the parameters used are listed in Table 0), H is given by 

1/2 

(9) 



H 



= 8.9-10" 



The disc is truncated if the viscous torque is smaller 
than the resonant torque. The criterion for the disc trun- 
cation at a given resonance radius is, in general, written 
as 



ml 



" ^CT m/ ) LR + ^C^mOcR < 0, (10) 
ml ml 



where the summation is taken over all combination of 
(m, I) which give the same resonance radius. For exam- 
ple, at the n:l resonance radius, at which Q — nOe, the 
summation is taken over all combinations of (to, I) with 
/ = n(rn — 1) for ILRs, I — n(m + 1) for OLRs, and I — nm 
for CRs. Actually, criterion (10) is determined only by the 



viscous torque and the torques from the ILRs of several 
lowest-order potential components, because the torques 
from the ILRs dominate those from the OLR and CR in 
circumstellar discs and high-order potential components 
contribute little to the total torque, even if the eccentric- 
ity of the orbit e is not small (Goldreich & Tremaine 1980; 
Artymowicz & Lubow 1994). 

Criterion ([To]) at a given resonance is met for a smaller 
than a critical value, a cr it- In Fig- @ 5 we plot a cr ;t at n : 1 
resonance radii. The resonant torques at the n : 1 radii are 
stronger than those at radii with other period commen- 
surabilities located nearby. It is important to note that 
«crit ^ 1 at the 2:1 and 3:1 resonance radii. The disc is 
truncated at the 3:1 radius if 0.17 ^ a <, 0.97 or at the 
4:1 radius if 0.037 <, a <, 0.17. Note also that the Be-star 
disc in 4U 0115+63 cannot be in a steady state. The gas 
supplied by the Be star will decrete outward and be ac- 
cumulated in the outermost part of the disc. This could 
contribute to make the disc dynamically unstable, as will 
be discussed in Paper II. 

It could be argued that the tidal truncation of the cir- 
cumstellar disc is an assumption implicit in our choice of a 
viscous decretion model for the disc, but it holds irrespec- 
tive of disc model as long as the outflow velocity in the disc 
is very subsonic, as observed for many Be stars. In order 
to show this, let us compare the timescale for disc trun- 
cation, Tt runc , with the drift timescale Tdrift- We can find 
Ttrunc via the viscous timescale t v [ s . Since the ratio of the 
viscous torque to the resonant torque is a/a cr ;t, we have 
nrunc ~ (a/a cri t)T vis - a~] t ( Ar/H ) 2 1 . In deriving the 
second similarity, we have used r v ; s ~ (Ar) 2 /ac s H ~ 
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Fig. 4. Critical values of a at n : 1 resonance radii. Anno- 
tated in the figure are the locations of the n : 1 commen- 
surabilities of disc and binary orbital periods. 



(Ar) 2 /aH 2 £l, where c s is the sound speed and Ar is a 
gap size (i.e., an interval between the disc outer radius 
and the radius at which the neutron star's gravity begins 
to dominate). On the other hand, the drift timescale is 
written as rarift ~ Ar/v r ~ M~ 1 (Ar/H)fl~ 1 , where v r 
and M. r are the radial velocity and Mach number, respec- 
tively. Consequently, we have 



Mr Ar 



Tdrift a C rit H 



(11) 



As shown in Fig. [|, a cr it ~ 1 at the 3:1 radius. From Eq.(g) 
with r/a ~ 0.5, we expect Ar/H = Ar/r/{H/r) ~ 0(1) 
for the disc of V635 Cas. Therefore, we conclude that the 
disc truncation timescale is much shorter than the drift 
timescale if the outflow is very subsonic. Note that this 
conclusion holds for Be/X-ray binaries in general, because 
all Be/X-ray binaries have similar disc parameters and 
the orbital parameters of these systems do not range very 
widely. 

4-4- Disc structure 

As mentioned above, the disc around V635 Cas cannot be 
steady. Nevertheless, it is instructive to study the steady 
disc structure as a model for the inner disc and/or the disc 
in the initial formation epoch. In this section, we assume 
the disc to be axisymmetric and azimuthally average the 
potential given by Eq. (fy). 



Okazaki (1997; see also Okazaki 2000b) showed that 
the equations which determine the velocity field in an 
isothermal decretion disc can be reduced to 



V r 



dV r 

dr 



GM* 



ad 



5cf_ 
2 r 



and 



£{R*) 



R> 



V r (R*) 



r 

V 



(12) 



(13) 



where £ = rV^ is the specific angular momentum, F ra( j 
is the vertically averaged radiative force, and V r and V$ 
are the radial and azimuthal components of the vertically 
averaged velocity, respectively. Eqs. ( |l2|) and ([l3]) show 
that the viscous decretion disc is an equatorial wind, in 
which material is slowly accelerated outward by the pres- 
sure force. 

Since the radial flow in Be-star discs is considered 
to be subsonic, the radiative force would arise not from 
the optically-thick strong lines but from an ensemble of 
optically-thin weak lines. Hence, for the radiative force, 
we adopt the parametric form proposed by Chen & Marl- 
borough (1994). 



GNU 



ad 



(14) 



where rj and e are parameters which characterize the force 
due to the ensemble of optically-thin lines. 

Given the radial velocity component at the stellar sur- 
face, the flow structure in the disc is obtained by solv- 
ing Eqs. ([l2]) and ( |l3| ) numerically. The surface density 
distribution is then obtained from the equation of mass 
conservation. In Fig. B, we show the structure of the 
isothermal decretion disc around V635 Cas with a = 0.1, 
V r (R*)/c s = 10~ 3 , and the outer radius at the 4:1 reso- 
nance. The parameters characterizing the radiative force, 
•q and e, are chosen so that the fundamental m — 1 mode 
in the disc has the period corresponding to the observed 
line profile variabilities (see Paper II). 

The solution shown in Fig. is not unique. We 
found that, for a — 0.1, there is a smooth solution for 
V r (R*)/c s < 3 x 10~ 3 . Any flow with a = 0.1 and 
V r (R*)/cs <; 3 x 10~ 3 , however, has a spiral critical point 
and does not go smoothly from the stellar surface to the 
disc outer radius. Thus, the radial flow in a decretion disc 
must be highly subsonic near the star. Such a flow satisfies 
the disc truncation condition, T trunc /rd r ift < 1, discussed 
in §[D| 

5. Discussion 

The truncation of the disc by the tidal/resonant interac- 
tion with the neutron star could explain many of the ob- 
served properties of V635 Cas. In N97, it was shown that 
large changes in the infrared magnitudes were not asso- 
ciated with significant changes in the associated colours. 
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Fig. 5. Structure of the viscous decretion disc around 
V635 Cas with a = 0.1 and K.(i?»)/c s = 1(T 3 . The 
radiative force in the form of rj(r/R*) e x GM^/r 2 with 
(77, e) = (0.1,0.1) is included. Solid, dashed, and dash- 
dotted lines denote a/a(R*), V r /c s , and V^/Vk(ii*), re- 
spectively, where Vk = (GM^/r) 1 / 2 is the Keplerian ve- 
locity. Two arrows denote the inner radius at r = i?» and 
the outer disc radius located at the 4:1 resonance radius. 

According to the models of Dougherty et al. (1994), this 
should correspond to a very small disc with a very high 
density - such that it is optically thick at all wavelengths 
for continuum radiation. The formation of such a disc is 
easy to explain if tidal truncation forces outflowing ma- 
terial to fall back on to the inner regions of the disc. It 
will also explain why the observed correlation between the 
strength of Ha and the circumstellar reddening observed 
in Be stars (Fabregat & Reglero 1990) does not seem to 
hold for Be/X-ray binaries - the disc could become very 
dense and the circumstellar reddening would not be mea- 
suring the size of the disc. Finally it would explain the cor- 
relation between the maximum strength of Ha and orbital 
period in Be/X-ray binaries found by Reig et al. (1997). 

The reality of this correlation has been questioned by 
Apparao (1998), who claims that the EW of Ha does not 
represent an appropriate measurement of the size of the 
disc and that the actual luminosity radiated in Ha does 
not show any correlation with the size of the orbit. We 
agree that the EW of Ha at a given time does certainly 
not represent a measurement of the size of the envelope. 
Many other factors such as the optical properties of the 
envelope, which seem to be highly variable, or an asym- 



-1 " 




Fig. 6. A description of the orbital model adopted for 
4U 0115+63, in the reference system centred on the Be 
star. The dash-dotted line represents the orbit of the neu- 
tron star. The thick dashed line represents the position 
of the first Lagrangian point (Li) around the orbit. The 
solid thin lines are the effective Roche lobes of the two 
stars at apastron and periastron (the position of the neu- 
tron star is marked with a cross). The labelled thick solid 
lines are the locations of the n : 1 commensurabilities of 
disc and binary orbital periods. Note that the radius of 
the 3:1 commensurability is smaller than the radius of the 
L\ point at periastron. 



metric density distribution in the disc (Negueruela et al. 
1998; see also Paper II) affect the observed EW. However, 
the maximum EW of Ha ever observed may provide a 
comparative estimate of the maximum size of an envelope 
if the source has been consistently monitored for a suffi- 
ciently long time - for example in the case of V635 Cas, 
where the typical timescales of variability are of the or- 
der of a few months (Paper II), ten years of monitoring 
provide a good estimate of the EW of Ha at maximum. 

Given that we expect the viscosity in the disc to be of 
the order of a ^0.1 as in accretion discs (e.g., Blondin 
2000; Matsumoto 1999), the disc around V635 Cas is likely 
to be truncated at the 4:1 resonance. Even if the viscosity 
is higher and becomes ~ 1, truncation still occurs at the 
3:1 resonance. In Fig. ^, it can be seen that even the 3:1 
resonance is significantly smaller than the radius of the 
first Lagrangian point (even when the neutron star is at 
periastron). Here, the potential ip describing the effects of 
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the gravitational and centrifugal forces on the motion of 
test particles orbiting the Be star is given by 



(15) 



where $ is the potential defined by Eq. (0). As a conse- 
quence, the neutron star will not be able to accrete a sig- 
nificant amount of material from the circumstellar disc. 
This fact explains the lack of series of Type I outbursts 
from 4U 0115+63. In order to have Type I outbursts in 
this system, we will need a large amplitude m = 1 mode 
to make the disc eccentric and cause Roche lobe overflow 
toward the neutron star. 

Accepting that the disc of Be/X-ray binaries is tidally 
truncated naturally explains why Be/X-ray binaries spend 
most of their time in quiescence. In general, the neutron 
star cannot accrete from the dense regions of the disc. The 
low-luminosity X-ray activity observed in systems such as 
A 0535+262 (Motch et al. 1991) would then be due to 
the accretion of some residual low-density outflow. In the 
case of V635 Cas, this material is not accreted because of 
the propeller mechanism. If we accept this picture, it is 
precisely the very existence of X-ray outbursts that needs 
to be explained. This will be addressed in Paper II. 

6. Conclusions 

By observing the source during a disc-loss episode, we have 
been able to determine the stellar parameters of V635 Cas. 
The object is classified as a B0.2Ve star at a distance 
of 7 — 8 kpc. Both the mass function and the estimated 
v sin i indicate a moderate inclination for the orbital and 
equatorial planes. The derived distance implies that the 
source can radiate close to the Eddington luminosity for 
a neutron star during bright outbursts. 

With the newly determined parameters, we have con- 
structed a model for 4U 0115+63. Based on the viscous 
decretion disc model for Be stars, we have numerically 
solved a criterion for tidal truncation and found that the 
disc surrounding V635 Cas must be truncated at a res- 
onance radius depending on the viscosity parameter and 
cannot be in a steady state. Although we have adopted 
a particular disc model, our conclusion is robust as long 
as the outflow in the disc is subsonic, a hypothesis sup- 
ported by several observational facts. Under normal condi- 
tions, the neutron star cannot accrete enough material to 
overcome the centrifugal barrier and switch on the X-ray 
emission. 
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